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Abstract 

New mode in the hadron scattering is predicted to appear at the 
energies beyond y/s — 2 TeV: the antishadow scattering mode and the 
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O i experiments at LHC and VLHC in hadronic reactions will be able to 

^ ■ reveal it. The appearance of the antishadow scattering mode at these 

\ energies is considered on the basis of unitarity and geometrical notions 

of hadron interactions. Connections with the nonperturbative-QCD 
models are discussed. 

X 

Introduction 

One of the most fundamental discoveries in hadron interactions at high ener- 
gies is the rise of total cross-sections with energy. It is accompanied by the 
rise of elastic and inelastic cross-sections as well as of the ratio of elastic to 
the total cross-section. 

For the first time the total cross-section increase was observed in K + p- 
interactions at the Serpukhov accelerator in 1970 [lj] and it was discovered 
later in pp-interactions at CERN ISR and at Fermilab || in other nucleon- 
and meson-proton interactions. Recent HERA data @] demonstrated the ris- 
ing behavior of the virtual photon - proton total cross-sections. Since then 
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a great progress in the experimental and theoretical studies of hadronic re- 
actions was achieved. Quantum Chromodynamics appeared as a theory of 
strong interactions and gave an explanation for the behavior of the observ- 
ables in the hard hadronic reactions, i.e. the reactions with high momentum 
transfers. However, the dynamics of long distance interactions (soft pro- 
cesses) is rather far from its understanding despite much work has been done 
in this field. The problems are directly related to the problems of confinement 
and chiral symmetry breaking. 

The approaches to soft hadronic processes are widely varied: Regge-type, 
geometrical or QCD-inspired models consider aspects of such processes from 
the different points of view and use various ideas on hadron structure and 
interaction dynamics. The major part of the models consider the global char- 
acteristics of hadron interactions such as a to t, Vinei and Udiff related to large 
distance interaction dynamics as reflecting gross features of hadron structure 
Despite of the difficulties in application of perturbative QCD for the 
description of long-distance interactions and their obvious nonperturbative 
character, it is often possible to represent the high-energy amplitude in the 
various model approaches as an expansion over a small parameter which de- 
pends on the kinematics of the process, e.g. for the case of non-increasing 
total cross-section the general form of the amplitude is 

F(s,t) = sJ2[r(s)] n exp 

n 

where r(s) ~ 1/ Ins is a small parameter at s — > oo. 

Since the expansion is not valid for the rising total cross-sections it is 
possible to find another representation for that case with the t-dependent 
expansion parameter |7J]: 

oo 

F(s, t) = sJ2 [T(V=£)] m $ m [R(s), MO, 

m=l 

where 




and Q m [R(s), yf—t\ is an oscillating function of transferred momentum. The 
above formulas as well as some other representations may be successfully 
used for the phenomenological analysis of the scattering amplitude at high 
energies. 



a(s)t 
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Thus, by now the theoretical treatment of soft hadronic reactions involves 
substantial piece of phenomenology and uses various model approaches. They 
are often based on divergent postulates, but their phenomenological parts are 
similar. In particular, an amplitude V(s,t) is considered as an input for the 
subsequent unitarization procedure: 

F(s,t) = $\V(s,t)]. 

To reproduce the total cross-section rise the input amplitude V(s,t) is usu- 
ally considered as a power function of energy. This function being taken as 
an amplitude itself violates unitarity in the direct channel. To obey unitarity 
in the direct channel an unitarization procedure should be used. 

There are several ways to restore unitarity of the scattering matrix. We 
consider two schemes: based on the use of eikonal and generalized reaction 
matrix respectively. There are also combined methods but those are not of- 
ten used. As it was mentioned various models for V(s, t) may be successfully 
used to provide phenomenological description of high energy hadron scatter- 
ing. However, in the particular model approaches the important dynamical 
aspects of interaction could be significantly obscured due to large number of 
free parameters. 

In this paper we discuss some general properties of hadron scattering, 
the implications of unitarity and analyticity, in particular, manifestations 
of the antishadow scattering mode and respective model predictions for the 
observables in elastic scattering and diffraction dissociation. Our main goal 
is to draw an attention to the existence of the antishadow scattering mode 
at the energies of LHC and VLHC. It might provide a new insight into the 
dynamics of diffraction and head-on hadronic collisions at superhigh energies. 

1 Geometrical Picture 

In the collisions of two high energy particles the de Broglie wavelength can be 
short compared to the typical hadronic size and hence optical concepts may 
be used as useful guidelines. Thus, the hadron scattering can be considered 
as a collision of two relativistically contracted objects of finite size. 

The relevant mathematical tool for description of high energy hadronic 
scattering is based on the impact parameter representation for the scattering 
amplitude. In the case of spinless particle scattering this representation has 
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the following form: 



g roo 

F(s,t) = -z bdbf(s,b)J (bV^i) 

7T Z JO 



(1) 



Note that for the scattering of particles with non-zero spin the impact pa- 
rameter representation for the helicity amplitudes has a similar form with 
substitution J — > J ax, where A A is the net helicity change between the final 
and initial states. The impact parameter representation as it was shown in 
|| is valid for all physical energies and scattering angles. This representation 
provides simple semiclassical picture of hadron scattering. 

It is often assumed, after the Chou-Yang model was proposed, that the 
driving mechanism of hadron scattering is due to overlapping of the two 
matter distributions of colliding hadrons. It could be understood by analogy 
with Glauber theory of nuclear interactions: one assumes that the matter 
density comes from the spatial distribution of hadron constituents and also 
assumes a zero-range interaction between those constituents. Such contact 
interaction might result from the effective QCD, e.g. based on the Nambu- 
Jona-Lasinio Lagrangian. 

The important role in the geometrical approach belongs to the notion of 
the interaction radius. The general definition of the interaction radius which 
is in agreement with the above geometrical picture was given in |J: 



where k = \fs/2 is the particle momentum in the c.m.s. The value for Iq(s) is 
chosen provided the contributions of the partial amplitudes from the angular 
momenta I > Iq(s) are vanishingly small. 

As a first approximation one can consider the energy independent inter- 
action intensity and describe the elastic scattering amplitude in terms of the 
black disk model where it has the form: 



Here R ~ 1/ is the interaction radius. The model is consistent with the 
observed structure in the differential cross-sections of pp- and pp-scattering 
at t close to 1 (GeV/cf. 

In the simplest case, neglecting the real part and spin, the impact param- 
eter amplitude f(s, b) can be obtained as an inverse transformation according 




(2) 
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to Eq. |I] with 

Thus, one can extract information on the geometrical properties of interac- 
tion from the experimental data. The analysis of the experimental data on 
high-energy diffractive scattering shows that the effective interaction area 
expands with energy and the interaction intensity — opacity — increases 
with energy at fixed impact parameter b. Such analysis used to be carried 
out every time as the new experimental data become available. For exam- 
ple analysis of the data at the ISR energies (the most precise data set on 
differential cross-section for wide t-range available for ^fs = 53 GeV) shows 
that one can observe a central impact parameter profile with a tail from the 
higher partial waves and some suppression (compared to gaussian) of low 
partial waves. The scattering picture at such energies is close to gray disk 
with smooth edge which is getting darker in its center with energy. 

Beside the above simple geometrical observations it is useful to keep in 
mind the rigorous bounds for the experimental observables. 



2 Bounds for observables and the experimen- 
tal data 

Bounds for the observables obtained on the firm ground of general principles 
such as unitarity and analyticity are very important for any phenomenological 
analysis of soft interactions. However, there are only few results obtained on 
the basis of the axiomatic field theory. 

First of all it is the Froissart-Martin bound that gives the upper limit for 
the total cross-section: 

a tot <C\n 2 s, (4) 

where C = 7r/m^ (= 60mb) and is the pion mass. 

Saturation of this bound, as it is suggested by the existing experimen- 
tal data, imply the dominance of long-distance dynamics. It also leads to 
number of important consequences for the other observables. For instance, 
unitarity leads to the following bound for elastic cross-section: 

/ \ \ a tot( S ) fr\ 
°el{s) > C 2 ■ (5) 

In s 
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Therefore, when the total cross-section asymptotically increases as In 2 s, 
elastic cross-section also must rise like In 2 s. It is important to note here that 
there is no similar bound for the inelastic cross-section and as we will see 
further the absence of such bound allows for appearance of the antishadow 
scattering mode at very high energies. 

If one considers a more general case when a tot oc In 7 s, then at asymptotic 
energies one should have 

ReF(s,0) ^ 7 vr 



ImF(s,0) 21ns 



and 



gjgtlgj a totK°J < ln -7/2( s ) (7) 

where of of (s) and cr" ot (s) are the total cross-sections of the processes a + b — > 
X and a+b — > X correspondingly. In the case of 7 = 2 the total cross-section 
difference of antiparticle and particle interactions should obey the following 
inequality 

Aa tot {s) < Ins. (8) 

Contrary to the total cross-section behavior, the existing experimental 
data seem to prefer decreasing Aa to t(s). Possible deviations from such be- 
havior could be expected on the basis of perturbative QCD fTof and it was 



one of the reasons for the recent discussions on the Pomeron counterpart - 
the Odderon. However, the recent measurements real to imaginary part ratio 
for forward pp scattering provide little support for the Odderon. We will not 
discuss more thoroughly the interesting problem of ReF/ImF ratio and will 
consider for simplicity the case of pure imaginary amplitude. 

For the slope of diffraction cone at t = in the case of a pure imaginary 
scattering amplitude the following inequality takes place: 



B(s) > -Mw (g) 
I8na el (s) 

which means that when the total cross-section increases as In 2 s, the same 
dependence is mandatory for the slope of diffraction cone. It is the stronger 
shrinkage than the Regge model predicts: B(s) ~ a' Ins. 



There is also bound [ 1 1 for the total cross-section of single diffractive 



processes. It was obtained by Pumplin in approach where inelastic diffraction 
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as well as elastic scattering are assumed to arise in form of a shadow of 
inelastic processes and has the form 



<Tdiff{s, b) < ^cr tot (s, b) - a ei (s, b). (10) 
The most significant assumption was that the diffractive eigenamplitudes in 



the Good-Walker []T2j picture do not exceed the black disk limit. 

At this point some details of the experimental situation have to be men- 
tioned. At the highest energies the experimental data for the total and elastic 
cross-sections, slope parameter of diffraction cone and cross-section of sin- 
gle inelastic diffraction dissociation have been obtained in pp-collisions at 
Fermilab. In particular, those measurements show that 

• the rise of the total cross-section of pp-interactions is consistent with 
In 2 s-dependence, however other dependencies are not ruled out; 

• elastic cross-section rises faster that the inelastic and total cross-sections 
and has a magnitude about 1/4 of the total cross-section. 

Comparing the value of the elastic to total cross-section ratio with the lower 
energy data one can conclude that the higher the energy, the higher both 
absolute and relative probabilities of elastic collisions. 

Impact parameter analysis [13|] of the data shows that the scattering 



amplitude is probably beyond the black disk limit \f(s,b) \ = 1/2 in head-on 
collisions. The Pumplin bound (Eq. [10]) is also violated in such collisions 
and this is not surprising if one recollects the original ad hoc assumption on 
the shadow scattering mode. 



3 Antishadow scattering mode 

The basic role in our consideration belongs to unitarity of the scattering 
matrix SS + = 1 which is a reformulation of the probability conservation. In 
the impact parameter representation the unitarity equation rewritten for the 
elastic scattering amplitude f{s,b) at high energies has the form 

Imf(s,b) = \f(s,b)\ 2 + r ] (s,b) (11) 
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where the inelastic overlap function r](s, b) is the sum of all inelastic channel 
contributions. It can be expressed as a sum of n-particle production cross- 
sections at the given impact parameter 

V (s,b) =5>n(*,&). (12) 

n 

As it was mentioned assumption of a pure imaginary amplitude is a rather 
common approximation at high energies and is adequate for our qualitative 



analysis. Then the unitarity Eq. 11 points out that the elastic scattering 



amplitude at given impact parameter value is determined by the inelastic 
processes. Eq. [11] imply the constraint 

LfM)l<i 

while the black disk limit presumes inequality 

|/M)|<l/2. 

The equality \f(s, b)\ = 1/2 corresponds to maximal absorption in the partial 
wave with angular momentum I ~ b^fsjl. 

The maximal absorption limit is chosen a priori in the eikonal method of 
unitarization when the scattering amplitude is written in the form: 

/(s,&) = i(l-exp[iw(s,&)]) (13) 

and imaginary eikonal u(s,b) = iQ(s,b) is considered. The function Q(s,b) 
is called opacity. Eikonal unitarization automatically satisfies the unitarity 
Eq. |ll] and in the case of pure imaginary eikonal leads to amplitude which 
is always obey the black disk limit. 

However, unitarity equation has the two solutions for the case of pure 
imaginary amplitude: 



f(s,b) = -[l±^l-Ms,b)]. (14) 

Eikonal unitarization with pure imaginary eikonal corresponds to the choice 
of the particular solution with sign minus. 

Several models have been proposed for the eikonal function. For instance, 
Regge-type models lead to the gaussian dependence of fl(s, b) on impact 
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parameter. To provide rising total cross-sections opacity should have a power 
dependence on energy 



n(s, b) oc s A exp[-b 2 /a{s)], (15) 

where a(s) ~ In s. In the framework of perturbative QCD-based models the 
driving contribution to the opacity is due to jet production in gluon-gluon 
interactions, when 

Q(s,b) oc cr Jet exp[— fib], (16) 

where <jj et ~ (s/s ) A . Such parameterizations lead to the rising total and 
elastic cross-sections and slope parameter: 

<?tot{s) ~ (?d{s) ~ B(s) ~ In 2 s (17) 

and the ratio 



O-el(s) 1 



(18) 



vtot(s) 2 

To include the mode where the scattering amplitude exceeds the black 
disk limit one should consider the eikonal functions with non-zero real parts. 
To ensure the transition from shadow to antishadow mode the real part of 
eikonal should gain an abrupt increase equal to tt at some s = Sq. The 
conventional models do not foresee such a critical behavior for real part of 
the eikonal. 

However, it does not mean that the eikonal model itself is in trouble. In 
particular, the account for fluctuations of the eikonal strongly modifies 



the structure of the amplitude and reduces it to algebraic form which is 
similar to that used in the unitarization scheme based on the generalized 
reaction matrix. 

The latter method is based on the relativistic generalization of the Heitler 
equation of radiation dumping fl5|j. In this approach the elastic scattering 



amplitude satisfies unitarity equation since it is constructed as a solution of 
the following equation [TjJ 

F = U + iUDF (19) 

presented here in the operator form. The Eq.|T^ allows one to satisfy unitarity 
provided the inequality 

Im£7(s, 6) > (20) 
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is fulfilled. The form of the amplitude in the impact parameter representation 
is the following: 

f( a b ) = U ^ b ^ (21) 
nS ' b) l-iU(s,b)' [ ' 

where U(s, b) is the generalized reaction matrix, which is considered as an 
input dynamical quantity similar to eikonal function. Similar form for the 
scattering amplitude was obtained by Feynman in his parton model of diffrac- 
tive scattering f]l|. Inelastic overlap function is connected with U(s,b) by 



the relation 

Construction of particular models in the framework of the [/-matrix ap- 
proach proceeds with the same steps as it does for the eikonal function, i.e. 
the basic dynamics as well as the notions on hadron structure are used to 
obtain a particular form for the [/-matrix. For example, the Regge-pole 
approach [[HJ provides the following form for the [/-matrix: 

U(s, b) oc is A exp[— b 2 /a(s)], a(s) ~ a']ns, (23) 

while the chiral quark model which will be discussed below gives the expo- 
nential 6-dependence 

U(s,b) oc zs A exp[-/ifo], (24) 

where \i is the constant proportional to the masses of the constituent quarks. 
We have mentioned here only the gross features of those model parameteri- 
zations without going into the details. 

The both parameterizations lead to In 2 s rise of the total and elastic cross- 
sections and slope parameter B(s): 

atot(s) ~ a el {s) ~ B{s) ~ In 2 a (25) 

at s — > oo. The above results are similar to conclusions of eikonal unitariza- 
tion. 

However, these two unitarization schemes lead to different predictions for 
the inelastic cross-sections and for the ratio of elastic to total cross-section. 
This ratio in the [/-matrix unitarization scheme reaches its maximal possible 
value at s — ► oo, i.e. 

- 1, (26) 
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which reflects in fact that the bound for the partial-wave amplitude in the 
[/-matrix approach is |/| < 1 while the bound for the case of imaginary 
eikonal is (black disk limit): |/| < 1/2. 

When the amplitude exceeds the black disk limit (in central collisions at 
high energies) then the scattering at such impact parameters turns out to be 
of an antishadow nature. It corresponds to the solution of unitarity equa- 
tion Eq. |ll] with plus sign. In this antishadow scattering mode the elastic 
amplitude increases with decrease of the inelastic channels contribution. 

The shadow scattering mode is considered usually as the only possible 
one. But the two solutions of the unitarity equation have an equal meaning 
and the antishadow scattering mode could also appear in central collisions 
first as the energy becomes higher. The both scattering modes are realized in 
a natural way in the [/-matrix approach despite the two modes are described 
by the two different solutions of unitarity Eq. [TJ|. 

Let us consider the transition to the antishadow scattering mode 



With conventional parameterizations of the [/-matrix in the form of Eq. £3| or 
Eq. |23| the inelastic overlap function increases with energies at modest values 
of s. It reaches its maximum value rj(s, b = 0) = 1/4 at some energy s = Sq 
and beyond this energy the antishadow scattering mode appears at small 
values of b. The region of energies and impact parameters corresponding to 
the antishadow scattering mode is determined by the conditions Imf(s, b) > 
1/2 and i](s,b) < 1/4. The quantitative analysis of the experimental data 
jn| gives the threshold value of energy: y/so ~ 2 TeV. 

Thus, the function rj(s, b) becomes peripheral when energy is increasing. 
At such energies the inelastic overlap function reaches its maximum value 
at b = R(s) where R(s) is the interaction radius. So, beyond the transition 
threshold there are two regions in impact parameter space: the central region 
of antishadow scattering at b < R(s) and the peripheral region of shadow 
scattering at b > R(s). At b = R(s) the maximal absorbtion (black ring) 
takes place (Fig. 1). 

The transition to the antishadow scattering at small impact parameters 
at high energies results also in a relatively slow rise of inelastic cross-section: 

^e,M=8,f | /_^ )|2 ~ln, (27) 

at s — > oo. 

It should be noted that appearance of the antishadow scattering mode 
does not contradict to the basic idea that the particle production is the driv- 
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b=R(s) b b=R(s) b 



Figure 1: Shadow and antishadow scattering regions 



ing force for elastic scattering. Indeed, the imaginary part of the generalized 
reaction matrix is the sum of inelastic channel contributions: 

ImU(s,b)=Y,U n (s,b), (28) 

n 

where n runs over all inelastic states and 

U n (s,b) = JdT n \U n (s,b,{U}\ 2 (29) 

and dT n is the n-particle element of the phase space volume. The functions 
U n (s, b, {£n}) are determined by the dynamics of 2 — > n processes. Thus, 
the quantity IraU (s, b) itself is a shadow of the inelastic processes. However, 
unitarity leads to self-damping of the inelastic channels J20j and increase of 



the function ImU(s, b) results in decrease of the inelastic overlap function 
T](s,b) when ImU(s,b) exceeds unity. 

At the energies when the antishadow mode starts to develop (it pre- 
sumably could already occur at the energies of the Tevatron-Collider) the 
Pumplin bound Eq. |10| for inelastic diffraction dissociation cannot be applied 
since the main assumption used under its derivation is not valid any more. 



4 The two modes of hadron scattering and 
the preasymptotic effects 

In this section we give a specific analysis of the hadron scattering on the 
basis of particular model. In Refs. [^l], ^] we the notions of effective chiral 
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quark model were used for the description of elastic scattering at small and 
large angles. Hadron dynamics is considered in the framework of effective 
Lagrangian approach. 

A common feature of the chiral models |HJ is the representation of a 
baryon as an inner core carrying the baryonic charge and an outer conden- 
sate surrounding this core [2l|. Following these observations it is natural to 



represent a hadron as consisting of the inner region where valence quarks 



are located and the outer region filled with quark condensate p2| . Such a 
picture for the hadron structure implies that overlapping and interaction of 
peripheral condensates in hadron collision occurs at the first stage. In the 
overlapping region the condensates interact and as a result virtual massive 
quarks appear. Being released part of hadron energy carried by the periph- 
eral condensates goes to a generation of massive quarks. Besides mass, quark 
acquires an internal structure and a finite size. Quark radii are determined 
by the radii of the clouds. Strong interaction radius of quark Q is determined 
by its Compton wavelength: 

r Q = Z/mQ, (30) 

where constant £ is universal for different flavors. In the model valence quarks 
located in the central part of a hadron are supposed to scatter in a quasi- 
independent way by the produced virtual massive quarks at given impact 
parameter and by the other valence quarks. 

The function U (s, h) (generalized reaction matrix) |T5| - the basic dy- 
namical quantity of this approach — is chosen as a product of the averaged 
quark amplitudes 

N 

U(s,b) = IlUeM)} (31) 



in accordance with assumed quasi-independent nature of valence quark scat- 
tering. The ^-dependence of the function (fg) related to the quark formfactor 
Fq{c[) has a simple form (/q) oc exp(— mg6/£). 

Thus, the generalized reaction matrix (in a pure imaginary case) gets the 
following form 



U{s,b) = ig 
where M = J2^=i m Q- 



/—\N 



ex P (-M6/0, (32) 
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At moderate energies s <C s (where ^/so = mg/a)the function U(s,b) 
can be represented in the form 



U(s,b) = ig 



1 + Na 



m Q 



exp(-Mfo/0- 



(33) 



At very high energies s > So ffe could neglect the energy independent term 
in (|32|) and rewrite the expression for U(s,b) as 

N/2 



U(s,b)=ig(s/m 2 Q ) exp(-M6/0. 



(34) 



Calculation of the scattering amplitude is based on the impact parameter 
representation and the analysis of singularities of F(s, (3) in complex /3-plane 

0- 

Besides the energy dependence of these observables we will emphasize its 
dependence on geometrical characteristics of non-perturbative quark inter- 
actions. 

The total cross-section has the following energy and quark mass depen- 
dencies 

< 2 



&tot{S) 



( m Q) 



Ms,N), 



(35) 



where (mg) = J2q=i m Q is the mean value of the constituent quark masses 
in the colliding hadrons. The function $ has the following behavior: 



<$>{s,N) 



ig/N 2 ) [1 + Na^~s/m Q } , s < s , 
In 2 s, s > s - 

Thus, at asymptotically high energies the model provides 

&tot(ab) 



(36) 



liU1 , 7 , 



1. 



Linear with yfs preasymptotic rise of the total cross-sections is in agree- 



ment with the experimental data up to ^fs ~ 0.5 TeV |19 



The inelastic cross-section can be calculated in the model explicitely, viz: 



Cinel(S) 



8ire 
N 2 (m Q y 



In 



1+0(1 + 



a-s/s. 
m Q ' 



N 



(37) 
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At asymptotically high energies the inelastic cross-section rise is as follows 

amM = W^y lns (38) 

At s ^> So the dependence of the hadron interaction radius R(s) and the 
ratio a e i/a to t on (itiq) is provided by the following equations: 

R(s) = ^b, (39) 

= (40) 
Ms) ^ ln s 

It is important to note here that such a behavior of the ratio o e il<Jtot and 



o'inei(s) results from self-damping of inelastic channels f2(J at small impact 
distances. Numerical estimates |19| show that the ratio u e i(s) / a to t(s) bo- 
comes close to the asymptotic value 1 at extremely high energies ^/s = 500 
TeV. 

Thus, unitarization drastically changes the scattering picture: at lower 
energies inelastic channels provide dominant contribution and scattering am- 
plitude has a shadow origin while at high energies elastic scattering domi- 
nates over inelasic contribution and the scattering picture corresponds to the 
antishadow mode. The functional s-dependencies of observables also differ 
significantly. For example, s-dependence of total cross-section at s <C so is 
described by a simple linear function of ^/s. It has been shown that such 
dependence does not contradict to the experimental data for hadron total 
cross-sections up to ^fs ~ 0.5 TeV. Such dependence corresponds to that of 
the hard Pomeron with A = 0.5, however, it was obtained in different ap- 



proach [f22| . This is a preasymptotic dependence and it has nothing to do with 
the true asymptotics of the total cross-sections. In the model such behav- 
ior of the hadronic cross-sections reflects the energy dependence of number 
of virtual quarks generated under condensate collisions in the intermediate 
transient stage of hadronic interaction. 

5 Antishadow scattering mode and inelastic 
diffractive processes 

Inelastic diffractive production as well as elastic scattering at low transferred 
momenta are the two basic processes which would lead to understanding of 
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large distance dynamics and hadron structure. Concerning inelastic diffrac- 
tive processes this statement can be traced back to the seminal paper |12| 



where such processes were considered as a result of a difference in absorption 
of various proton states. Later on these states have got a parton-like inter- 
pretation. New data were obtained for single diffraction production process 

h 1 + h 2 -*h 1 + h* 2 (41) 

when the hadron h 2 is excited to the state h 2 with invariant mass M and 
the same quantum numbers. Its subsequent decay results in the multiparticle 
final state. The inclusive differential cross-section shows a simple dependence 
on the invariant mass 



d(Tdiff 



« TFf (42) 



dM 2 M 2 ' 

However, energy dependence of the diffractive production cross-section <7diff( s ) 
is not so evident from the data. This ambiguity is partly due to difficulties 
in the experimental definition of the inelastic diffractive cross-section. 

The particular experimental regularities observed in diffractive produc- 
tion can be described in the framework of different approaches. 1/M 2 depen- 
dence is naturally described by the triple-pomeron diagrams in the frame- 
work of Regge-model. The proposed in Ref. [25| similarity between the 



Pomeron and photon exchanges allowed to calculate diffractive dissociation 
cross-section in terms of structure function vW 2 measured in deep inelastic 
lepton scattering. Several models use optical picture for the description of 
diffractive production |26| but these models in large extent concern the angu- 
lar distribution of diffractive cross-section and M 2 -dependence is left beyond 
of their scope. The attempt to explain M 2 -dependence in the framework of 
optical model considering diffractive dissociation as a bremsstrahlung where 
virtual quanta are released from a strong field was made in Ref. |27| . 

In this section for description of single diffractive processes we use model 
approach described in section 4. 

To obtain the cross-section of the diffractive dissociation process we should 
single out among the final states in Eq. those corresponding to the process 
(|4l| ) . Let for simplicity consider again the case of pure imaginary [/-matrix. 
Then we can represent da,n ff/ dM 2 in the following form 

-dNp-^h Mb [l + U(s,b)] 2 { ] 
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where expression for L^i//(s, b, M) includes contributions from all the final 
states \n)diff which results from the decay of the excited hadron h 2 of mass 
M: h* 2 -> \n) dif f. 

For consideration of the diffractive production at the quark level we ex- 
tend the picture for hadron interaction for elastic scattering, described in 
section 4. Since the constituent quark is an extended object there is a non- 
zero probability of its excitation at the first stage of hadron collision during 
the interaction of peripheral condensates. Therefore it is natural to assume 
that the origin of diffractive production process is the excitation of one of 
the valence quarks in colliding hadron: Q — > Q*, its subsequent scattering 
and decay into the final state. The excited constituent quark is scattered 
similar to other valence quarks in a quasi-independent way. The function 
Udiff(s,b,M) can be represented then product 

N-l 

U diff (s,b,M) = (f Q *(s,b,M Q *)) J] (/qM)>, (44) 

0=1 

where Mq* is the mass of excited constituent quark, which is proportional 
to the mass M of excited hadron h 2 for large values of M. The last state- 
ment presumes the additivity of constituent quark masses. The 6-dependence 
of the amplitude (/q*) is related to the formfactor of excited quark whose 
radius is detemined by its mass Mq* [tq* = £/Mq*). The expression for 
Udiff{s, b, M) can be rewritten then in the following form: 

U&ffis, 6, M) = g*U(s, b) exp[-(M Q * - m Q )6/£], (45) 

where constant g* is proportional to the relative probability of excitation of 
the constituent quark. The value of g* is a non-zero one, however, g* < 1 since 
we expect that the excitation of any constituent quark has lower probability 
compared to probability for this quark to stay unexcited. The excited quark 
is not stable and its subsequent decay is associated with the decay of excited 
hadron h* 2 into the multiparticle final state \n)diff- 

The cross-section of diffractive dissociation process is given by expression 
( |43| ) and has the following s and M 2 dependence 

llAP " (M Q .-m^ (a ' 0) " -M^^ 0) (46) 

Thus, we obtained the familiar 1/M 2 dependence of the diffraction cross- 
section which is related in this model to the geometrical size of excited con- 
stituent quark. 
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The double dissociation processes 



h 1 + h 2 ^h* 1 + h* 2 (47) 

can be considered on the grounds of previous approach to the single diffractive 
dissociation. Here one of the constituent quarks in each of the colliding 
hadrons should be excited. Cross-section of double diffraction process has 
similar M 2 - and s-dependencies and is to be suppresed in comparison with 
the single diffractive cross-section by an extra factor g* < 1. 

The energy dependence of single diffractive cross-section has the following 
form 

* 9 / f M ? dM 2 M ~, ., sfl-iti) 

a diff (s) = 8ng*ev(s,0) J M? -j^r = ferg'^fo °) ln M< 2 , ( 4 §) 

where X\ is the lower limit of the relative momentum of hadron h\{x\ ~ 0.8 — 
0.9) which corresponds to the experimental constraint on diffractive process. 
Eq. ( Ji8"D shows that the total cross-section of diffractive dissociation has a 
non-trivial energy dependence which is determined by the contribution of 
inelastic channels into unitarity equation at zero value of impact parameter. 
The dependence of 77 (s, 0) is determined by Eq. (0), where expression for 
U(s,b) is given by Eq. fl52]). At s < s , (s is determined by equation 
|C/(so)0)| = 1) T](s, 0) increases with energy. This increase as it follows from 
Eq. (|32|) and from the experimental data |28| is rather slow one. However at 
s > s , 77(s,0) reaches its maximum value 77 (s, 0) = 1/4 and at s > s , the 
function r](s, 0) decreases with energy. At s —> 00: 

( 1 V 

(7 di ff(s) oc I -j= I Ins (49) 

since 77(5,0) oc (l/y/s) N in this limit. 

Thus at asymptotical energies the inelastic diffraction cross section drops 
to zero. Decrease of diffractive production cross-section at high energies 
(s > s ) is due to the fact that 77(5, b) becomes peripheral at s > s and the 
whole picture corresponds to the antishadow scattering at b < R(s) and to 
the shadow scattering at b > R(s) where R(s) is the interaction radius. The 
qulitative behavior of cr<nff(s) is shown on Fig. 2. 

The development of the antishadow mode in head-on pp and pp-collisions 
could be associated with new phenomena in the central hadronic collisions 
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Figure 2: Energy dependence of diffractive cross-section 

where the temperatures are high and the energy density can be up to several 
GeV/ fm 3 . In such collisions the constituent quarks have noticeable prob- 
ability to be excited. Due to its high mass and small transverse size the 
excited state has low probability of interactions with other particles. It may 
be also related to an interesting phenomena in cosmic ray experiments where 
particles with abnormal persistency in lead chambers were observed p9[| . 

Of course, there might be different reasons for decrease of adiff(s). The 
decreasing energy dependence of (Jdiff{ s ) was also predicted in Refs. |30| , 



3jJ . As it was pointed out in Ref. [|12| in the limit of complete absorption 



the diffractive dissociation should vanish. It was advocated in Ref. |TJj 
that this situation will occur at superhigh energies and it is the reason for 
decrease of inelastic diffractive cross-section. This is completely the same 
behavior as it is predicted by the model presented, however in our case the 
reason for that is the transition to the antishadow scattering mode in head-on 
collisions in the multi-TeV energy range. It should be noted, however, that 
the diffractive cross-section at preasymptotic energies has a similar to total 
and elastic cross-section energy dependence and it will be discussed in the 
concluding part of this paper. 

6 Universal preasymptotics 

The straitforward interpretation of the recent HERA data on the deep- 
inelastic scattering together with the analysis of the data on hadron-hadron 
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scattering in terms of the Regge model could lead to the unexpected conclu- 
sion on the existence of the various Pomerons [[32|] or the various manifesta- 



tions of unique Pomeron in the different processes depending on the typical 



scale of the process |33] . The approaches |34|, |35| contending the dominance 
of the soft Pomeron do not rule out existence of the hard Pomeron either. 
Indeed, soft hadronic reactions imply that the Pomeron's intercept a-p = 



1.08 ||32| , and small-x dependence of the structure function F 2 (x,Q 1 ) leads 



to a-p = 1.4-1.5 |36| [37] and the measurements of the diffractive cross-section 



in the deep inelastic scattering provide a-p = 1.23 [p8| . So, does this mean 
that we have few Pomerons or we have few different manifestations of the 
same Pomeron depending on the particular process? Probably both options 
are not to be considered as the firm ones, since the experimental data used 
to advocate these statements were obtained at not high enough energies 
where, in fact, the preasymptotic regime of interactions does take place. The 
above conclusions are based on the presumed dominance of the Pomeron 
contribution already in the preasymptotic energy region. What is called 
a Pomeron is to be interpreted as a true asymptotical contribution of the 
driving mechanism. 

In this section we argue that all the three classes of the processes de- 
scribed above are related to the similar mechanisms and the corresponding 
energy dependence of the cross-sections can be well described by the universal 
functional energy dependence of the type a + byfs. Such dependence is valid 
for the preasymptotic energy region only and beyond this region unitarity 
changes the picture drastically. We consider for illustration the unitarized 
chiral quark model (section 4). 

Fit to the total hp cross-sections gives small values for the parameters 
g and a ( g, a <C 1) fll9|. It means that at s <C Sq the second term in the 



square brackets in Eqs. (^TJ) and (|2"2"| ) is small and we can expand over it. The 



numerical value of s is determined by the equation \U(s, 0)| = 1 and is |19 
~ 2TeV. At this energy the amplitude has the value |/(s ,0)| = 1/2. 
The value of So is on the verge of the preasymptotic energy region, i.e. the 
Tevatron energy is at the beginning of the road to the asymptotics. Evi- 
dently the HERA energy range W{= y/s^) < 300 GeV is in a preasymptotic 
domain. 

The above model gives the linear with yfs dependence for the total cross- 
sections according to Eqs. (|2l|) and (^2|): 

o&» = a + b v Ts, (50) 
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where parameters a and b are different for different processes and the same 
is true for the scale s . It was shown |Tj| that Eq. |5(] is in a good agreement 



with the experimental data. 

The same dependence for the total cross-section of 7*p scattering is as- 
sumed by the small-x behavior of the structure function F 2 (x,Q 2 ) |36|, |37 
and obtained in [3S]: 

F 2 (x, Q 2 ) = a(Q 2 ) + b(Q 2 )/^. (51) 

The experimental data also indicate the critical behavior of the function 
b{Q 2 ) at Q 2 ~ 1 (GeV/c) 2 . This scale could be related to the radius of a 
constituent quark and its structure. 

The third value for the Pomeron intercept a-p = 1.23 has been obtained 
from the analysis of the experimental data on the diffractive cross-section in 
deep-inelastic scattering |38|] where the dependence of do d ^^ XN j dM\ on W 
was parametrized according to the Regge model and the Pomeron dominance 
has been assumed: 

dcr%L XN /dM 2 x oc (W 2 ) 2a -- 2 . (52) 

The data demontrate linear rise of the differential cross-section da^l p {^ XN /dM x 
with W, i.e. we observe here just the same functional dependence on the 
c.m.s. energy as for o"^' 7P ' 7 p . Regarding the preasymptotic nature of the 
interaction mode we arrive to the universal c.m.s. energy dependence in the 
framework of the used model. 

Indeed, in the framework of this model the hadron inelastic diffractive 



cross-section is given by the following expression |4C 



dM\ ~ Ml V[S ^ h {b6} 

where 

rj(s, b) = lmU(s, b)/\l - iU(s, b)} 2 

is the inelastic overlap function. 

At the preasymptotic energies s sq the energy dependence of inelastic 
diffractive cross-section resulting from Eq. ( p2[) is again determined by the 
generic form 

oc a + ovs. (54) 



dM 



2 
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Inelastic diffractive cross-section for the 7*p interactions can be obtained 
using for example VMD model, i.e. 



datLL 



7 P- 



dM\ 



oc a{Q 2 ) + b(Q 2 )W. (55) 



The same functional dependence can be obtained using the "aligned jet" 



model |4T] along with the unitarized chiral quark model 

The above linear dependences for the cross-sections of different processes 
is the generic feature associated with the preasymptotic nature of the inter- 
action dynamics at s <C sq. As one goes above this energy range the function 
\U(s, 6)| is rising and when \U(s, 0)| > 1 the unitarity starts to play the ma- 
jor role and provides the In 2 s rise of the total cross-sections at s ^> s 
and also the following behavior of the structure function F 2 (x, Q 2 ) 

F 2 {x,Q 2 ) ocln 2 (l/x) (56) 



at x — > p9f . At the same time unitarity leads to the decreasing dependence 



of the inelastic diffractive cross-section at s — > oo 

(57) 




dM 



for the hp, 7p and Yp processes [0. Eq. |57| is associated with the anti- 
shadow scattering mode which develops at small impact parameters at s > sq. 

Thus, we might expect the different asymptotic and universal preasymp- 
totic behaviors for the different classes of the diffraction processes. 

To summarize, we would like to emphasize that the unified description 
of the processes of hp, 'jp and 7*p diffraction scattering with the universal 
cross-section dependence on the c.m.s. interaction energy is possible. For 
the illustration we used the unitarized chiral quark model which has a non- 
perturbative origin and leads to the linear c.m.s. energy dependence of the 
cross-sections in the preasymptotic energy region for the above processes. 
Universality of such preasymptotic behavior agrees with the experiment. 

The assumption on the existence of the different Pomerons results from 
the use of the asymptotic formulas in the preasymptotic energy region and 
the neglect of the unitarity at higher energies beyond this preasymptotic 
region. It should be taken with certain caution. 
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Conclusion 



Studies of soft interactions at the highest energies can lead to the discover- 
ies of fundamental importance. The genesis of hadron scattering with rising 
energy can be described as transition from the grey to black disk and even- 
tually to black ring with the antishadow scattering mode in the center. Such 
transitions are under control of unitarity of the scattering matrix. 

The appearance of antishadow scattering mode could be revealed per- 
forming impact parameter analysis of elastic scattering and directly in the 
measurements of the inelastic diffractive cross section (cf. Figs. 1,2). 

It would be interesting to speculate on the particular physical origin of 
the antishadow scattering mode. Its existence can be correlated with the new 
phenomena expected at high energies in the central hadronic collisions. Such 
collisions are usually associated with the formation of quark-gluon plasma 
and disoriented chiral condensate in the inner part of the interaction region. 
What are the particular correlations between those phenomena and the anti- 
shadow scattering? The answer can be obtained in the nonperturbative QCD 
studies and in the experiments devoted to studies of soft processes at LHC 
and VLHC It seems that the anomalies observed in cosmic ray experiments 
211 might also be correlated with development of the antishadow scattering 



mode in the central hadron collisions. 
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